Local exchange-correlation vector potential with memory in Time-Dependent Density 
Functional Theory: the generalized hydrodynamics approach 
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Using Landau Fermi liquid theory we derive a nonlinear non-adiabatic approximation for the 
exchange-correlation (xc) vector potential defined by the xc stress tensor. The stress tensor is 
a local nonlinear functional of two basic variables - the displacement vector and the second-rank 
tensor which describes the evolution of momentum in a local frame moving with Eulerian velocity. 
For irrotational motion and equilibrium initial state the dependence on the tensor variable reduces 
to that on the metrics generated by dynamical deformation of the system. 
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Density-Functional Theory (DPT) offers an alterna- 
tive to the common wave function or Green's functions 
methods for studying the properties of inhoniogeneous 
many-body systems. Whereas the equilibrium DFT has 
become a standard computational tool in solid state, 
atomic, molecular and nuclear physics^, the Time De- 
pendent DFT (TDDFT)2 is still intensively developing. 
An explosion of interest in TDDFT is not surprising since 
the traditional many-body methods encounter enormous 
computational difficulties if one attempts to consider re- 
alistic non-equilibrium systems. The central problem of 
any DFT is to find a good approximation for the xc po- 
tential. In the static DFT the construction of approxima- 
tions usually starts from the Local Density Approxima- 
tion (LDA), which, by itself, often gives very good results. 
A similar approximation in TDDFT is still lacking. The 
commonly used Adiabatic LDA (ALDA) can be justified 
only for systems with the energy gap A if the inverse 
characteristic time of the process = w ^ A. This 
condition is strongly violated in most practical cases, 
which makes a construction of the non-adiabatic func- 
tionals highly desirable. 

In the mid-nineties it was realized that any consistent 
incorporation of non-adiabaticity in the xc potential in- 
evitably destroys its spatial locality'^. An important step 
was made by Vignale and Kohn (VK)'* who showed that 
in the linear response regime switching of the basic vari- 
able from density n to current j allows to construct a 
consistent local non-adiabatic approximation. The VK 
approximation is exact for Lkp ^ 1 (which is a standard 
quasiclassical LDA condition, where L is the length scale 
of the density variation and kp \s the local Fermi mo- 
mentum), and ujL/vf 2> 1. Vignale, Ullrich and Conti 
(VUC)^ realized that velocity v = j/n is a more nat- 
ural variable than the current, as it allows to cast the 
VK result in a transparent visco-elastic form. Yet avail- 
able up to now nonlinear non-adiabatic functionals^''' 
have been designed phenomenologically taking into ac- 
count a number of constraints due to the harmonic po- 
tential theorem^, zero force'^ and torque condition and 
the known VK linear form. However, all these require- 
ments, which must certainly be met, do not uniquely 



determine the nonlinear xc functional. 

The non-adiabatic LDA is in fact closely related to the 
generalized hydrodynamics, which we derived in Ref. 8 to 
remove some deficiencies of the adiabatic Bloch's theory^. 
In this paper we extend the hydrodynamics formalism to 
the full-interacting case and derive the nonlinear non- 
adiabatic local approximation for the xc vector poten- 
tial. The hydrodynamical approximation becomes exact 
under the VK conditions Lkp ^ 1, LoLjvp ^ 1, plus 
the additional requirement that ep j lo \ {ep \% the lo- 
cal Fermi energy), which is a temporal analogue of the 
inequality Lkp ^1. In the linear regime our local xc 
functional transforms to the VUG form. It also shows 
an initial state dependence, inevitable in any non-linear 
TDDFT^". 

Let us first formulate TDDFT in the form which 
is convenient for further discussion. Gonsider a non- 
equilibrium many-body system with Hamiltonian H = 
T -I- t/cxt + W^, where T is the kinetic energy operator, 
C4xt — J djdJgxt{'^,t)n{'x.) describes the interaction with 
the external potential, and W is the interaction Hamil- 
tonian related to the interparticle interaction Vx.-x' ■ 

The averaged equations of motion for density n and 
current j operators can be represented in the form of the 
hydrodynamical conservation laws 



Dtn 
mnDtv^ - 



"-<9pVp = 0, 



(1) 
(2) 



where Dt — dt + v^di, is the co- moving derivative and 
U = Ucxt+Un is a sum of external and Hartree potentials. 
Since the Hartree term is singled out in Eq. (2), the stress 
tensor P^i, — T^^ + W^i, contains only the kinetic T^^ and 
the xc W^^ parts^^ 
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T,u = ^((7r;^+)(7r.^) + (7r:^+)(^^^) - ^V^n)(3) 



+Ax';x-(1-A)x 

,d\, (4) 



where tt^ — 



-idfj_ — mv^ is the relative momentum. 



Gx;x' — (V''''(x)n(x')'0(x)) — n(x)n(x') is the pair cor- 
relation function and the angular brackets stand for ex- 
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act many-body state average. We use the summation 
convention over repeated indeces throughout the paper. 

According to the first part of the Runge-Gross (RG) 
proof of TDDFT^^ (see also^'^), the many-body wave 
function \E'(f) with initial condition \l/(0) =^"05 is a func- 
tional of the current density j or, equivalently, of the 
velocity v. Hence the stress tensor is a functional of the 
velocity and of the initial many-body state P^i/[^0)V]. 
Therefore Eqs. (1) and (2) are formally closed and RG 
theorem can be viewed as a proof of existence of the exact 
quantum hydrodynamics with the memory of the initial 
many-body correlations'^^. 

The important part of DFT is the Kohn-Sham (KS) 
construction, which, in a dynamical case, can be intro- 
duced as follows. Let us consider a system of noninter- 
acting KS particles moving in the self-consistent scalar 
U = C/oxt + Uh and vector A^'^ potentials. The local 
conservation laws for this system take the form 

Dtn + nd^v^=0, (5) 
mnDtv^ + d^P^^-nF^'' + nd^U = 0, (6) 

where F'^'^ = — 9i A'^^+vx (Vx A'^^) is the force due to the 
xc vector potential, and KS stress tensor is defined 
by Eq. (3) with tt^ = — — mw^ — A^. Comparing 
Eqs. (5), (6) with Eqs. (1), (2) we find that the density n 
and the velocity v of KS and interacting systems coincide 
if the functional A^'^[n,v] satisfies the equation 

F^' = -f (v X (V X A-))^ = -^5.AP-, (7) 

(8) 



A p^^ = p — P^ 



Equation (7) is exactly of the form suggested by VUC^ 
to satisfy Newton's third law'^. Apparently A'^'^ is a 
functional of the interacting initial state \E'o (via P^j/) and 
the KS initial state ^'q (via P^^) (for a recent discussion 
of the initial state dependence see Ref. 10). 

Let us now turn to the local approximation for A^'^. 
Equation (7) reduces the problem to the construction of 
approximation for the xc stress tensor AP^^. As usual 
in LDA, we require that Lkp ^> 1. In dynamics it is 
natural to assume the time-analogue of this condition: 
3> 1. At low temperatures the two conditions bring 
us to the domain of applicability of Landau Fermi-liquid 
theory'^. The latter is defined by the energy functional 
E\np] (local in space and time) of the quasiparticle dis- 
tribution hmction n-p (x, t) . An explicit example of this 
(in general nonlinear) functional is given by the x-only 
approximation 



(9) 



The function rip satisfies the Boltzmann-type kinetic 
equation^^. It is convenient to separate the convective 

motion with the velocity v and to introduce the distribu- 



tion function of relative momentum n-c 



^p-i-mv? 



which 



satisfies the following equation^ 



9ep drip 



mDtVv + 



dU 



drip 



/[np] (10) 



with initial condition np(x, 0) = A^(p, x). Here /[rip] is a 
collision integral, and the quasiparticle energy is defined 
as Ep = SE/Srip. Equation (10) differs from the com- 
mon kinetic equation'"'' in two points. First, the time 
derivative dt is replaced by the co-moving derivative Dt 
and, second, the force term contains additional inertial 
contributions (the first two terms in the square brack- 
ets). The latter arise due to the transformation to the 
local frame that moves with Eulerian velocity v. The 
first term mDtVu is the linear acceleration force whereas 
the second term p^d^Vi, is a sum of a force due to the 
deformation rate \p^{dvVn + dnVv) and the Coriolis force 
\Pti{dtiV^ - d,yVf,). 

Equation (10) leads to the local conservation laws of 
Eqs. (1), (2) with the stress tensor expressed in terms of 
the Landau functional'^ 



p p 

In the x-only approximation Eq. (11) reads 



m 



p>p' 



dpv 



+ 



2 P- 



(11) 



UpTlpi 



(12) 



Substituting d^U from Eq. (2) into Eq. (10) we eliminate 
the explicit dependence on the external potential from 
the kinetic equation: 



Dtfi 



dpv dxv 



dxn 



dP, 



der. 



dxv ndxu 



dn 
dpu 



P =7 



(13) 

Equations (13) and (11) taken together with the defini- 
tion of the density n = hp constitute a closed uni- 
versal (i.e. independent of external potential) problem. 
For a given velocity v(x, t) and initial conditions, this 
problem defines a universal functional P^i^[iV(p, x), v]. 

Calculation of the xc stress tensor AP™ of Eq. (8) 
requires a knowledge of the KS stress tensor P^^. The 
latter is obtained by the repetition of the above consid- 
eration for the KS system. The resulting noninteracting 
universal problem is 



Dtni 



Pu 



m dxv 



dvi, 
dx,. 



(p X (V X A-)), 



ndxu 



dhl 

OPu 



(14) 



where rip is a distribution function of KS particles and 
the KS stress tensor is defined as follows 



E 



/(.p. 



(15) 
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Equations Eqs. (13), (11), (14), (15), and Eqs. (7), (8) 
allow to construct the dynamical counterpart of LDA, 
provided the Landau functional -E[np] is known. 

Since the exchange self-energy is local in time, in x-only 
case the collision integral in Eq. (13) vanishes, and we can 
ignore the condition £f/w 3> 1. Thus, irrespectively of 
the frequency, A'^ obtained with the Landau functional 
of Eq. (9) is exact if Lkp ^ 1- This is a direct dynamical 
analogue of a static x-only LDA potential v^-^"^, which 
provides the long wave length limit of the time-dependent 
Optimized Effective Potential (OEP) approximation. 

In general, Eqs. (13), (11) seem to be intractable. They 
can be, however, solved in some important limiting cases. 
In a collision-dominated regime Vc/ J^3ix{u),vf/L} » 1 
{vc is a collision frequency) the memory is rapidly lost 
and the solution recovers the classical Navier- Stokes form 
of P^i^^^. On the contrary, the memory effects are pro- 
nounced in the opposite limit of a fast coUisionless mo- 
tion. In this regime the solution to Eqs. (13), (11) can 
be found analytically if Lcu/vp 3> 1, which is the VK 
condition^ and, at the same time, the condition of appli- 
cability of the generalized hydrodynamics^. In the rest 
of the paper wc analyze this case in detail and derive the 
local non-adiabatic approximation for AP^^. 

Let us estimate different terms in Eq. (13). The time- 
derivative term and the first term in brackets in Eq. (13) 
are of the order of ui (the continuity equation requires 
v/L ^ Lu), whereas all other terms give contributions 
~ Vp/L. Hence to the leading order in vf/ijjL <C 1 only 
inertial forces are relevant, and, therefore, the function 
np satisfies the following equation 



d 



d 



dvu d 



dt ^'^ dx„ dxu dp„ 



hp{-x.,t) = 



(16) 



with initial condition np(x, 0) = iV(p,x). In Eq. (16) we 
used an explicit expression for the co-moving derivative 
Dt. With the same accuracy, the KS distribution func- 
tion np(x, t), which defines the KS stress tensor of 
Eq. (15), also satisfies Eq. (16), but with another initial 
condition np(x, 0) = N^{p,x). 

To solve Eq. (16) we introduce a nonlinear transforma- 
tion of variables x, p ^ ^, k: 



x = x(ei) =C + u($,f), 
P = p(^,k,t), 



(17) 
(18) 



where the functions x(^, t) and p(^, k, t) are solutions to 
the following initial value problems 



ax(t) 

dt 



:v(x(t),i), X(0)=C 



(19) 



«)=-^!VMMp,(,), ,^(o) = fc,. (20) 

In terms of the new variables Eq. (16) transforms to 
the equation 9tn(k, ^,t) = 0, which is trivially solvable. 
Hence the solution of the original problem is 



where ^(x, t) and k(x, p,t) are obtained by inversion of 
Eqs. (17) and (18). Variable ^, Eq. (17), is a Lagrangian 
coordinated^ which has a meaning of the initial position 
of the fluid element presently at x. Vector u = x — ^ 
is thus a displacement of the fluid element. Similarly, k 
is the initial momentum of a quasiparticle initially at ^ 
that moves under inertial forces and, at the time moment 
t, acquires a momentum p. 

In the Eulerian (spatial) description^^ the above solu- 
tion for the quasiparticle np(x, t) and KS np(x, t) distri- 
bution functions can be written as 

np(x, t) = iV(p,?7;J(x, t), X - u(x, t)), (22) 
n^\x, t) = iV^(p„r7;J(x, t), x - u(x, t)), (23) 

where the displacement vector u(x, t) and the second- 
rank tensor 77^y(x, t) are defined as follows 

Au(x,i)=v(x,t), u(x,0)=0, (24) 

A?y/*i/(x,i) = -^^r?a,.(x,t), r]u.„{x,0) = S^^ (25) 

Tensor rj^^ relates the momenta p and k and, therefore, 
describes the rotation and the stretching of a quasiparti- 
cle momentum, which are caused by the inertial forces in 
the frame moving with Eulerian velocity. Using Eq. (25) 
one can prove that this tensor has an important prop- 
erty - the determinant of the matrix t]^i, equals to the 
Jacobian J of the coordinate transformation of Eq. (17). 
Hence detrj^i, is directly related to the density n(x, t): 



det V,. = det 1^ = J(x, t) = ^^^4^, 
ox^ no($(x,i)) 



(26) 



where r?.o(x) is the initial density distribution. In a par- 
ticular case; of irrotational motion the Coriolis force van- 
ishes (5^ IV = SiyUp), and the solution to Eq. (25) is lo- 
cally expressed in terms of deformation gradients 



dxu, 



dXa 



(27) 



np(x,i) =iV(k(x,p,<),^(x,i)), 



(21) 



Substituting the distribution functions Eqs. (22) and (23) 
into the definitions of the stress tensors P^^, Eq. (11), and 
P^^, Eq. (15), and calculating their difference, Eq. (8), 
we obtain the final expression for AP^^. The xc stress 
tensor APJ^ is a local functional (function) of two basic 
variables ^(x, f) = x — u(x, t) and rj^vix, t) 

AP-(x, t)) = l^P^lMx, i), ^(x, t)). (28) 

It is also a functional of the initial distribution functions 
iV(p, x) and N^{p, x) as it must be for any consistent ap- 
proximation in TDDFT. If the evolution starts from the 
ground state we have N = TV''' = fp{x) = 6{kp{x) — p'^), 
and the initial state dependence reduces to the depen- 
dence on the initial density no(x) (here kpix) is the lo- 
cal Fermi momentum at t = 0) . In addition, the isotropy 
in p-space of the equilibrium initial distribution function 
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requires that AP™ should not depend directly on 77^,^, 
but only via the symmetric tensor g^^, = rj^aVi^a' 

AP-(x,t)) = AP-(5„^(x,i),no(C(x,i))). (29) 

Interestingly, in a special case of irrotational motion, 
using Eq. (27) we obtain g^^{x,t) = ^Hj, which is 
exactly the metrics generated by the coordinate trans- 
formation of Eq. (17): (rf^)^ = g^vdx^dxi,. 

In the linear response regime Eq. (22) leads to the fol- 
lowing expression for the deviation of the distribution 
function rip from its initial value iV(p, x) 

Substituting Eq. (30) into Eq. (11) and considering the 
equilibrium initial distribution iV(p,x) = /p (x), we re- 
cover the correct elastic form of the stress tensor^. We 
note, that viscosity vanishes within the limits of applica- 
bility of Landau theory. 

The expression for the xc stress tensor Eq. (29) obvi- 
ously satisfies all general theorems and limiting require- 
ments, as initial microscopic equations do. It is note- 
worthy that the dependence on a delayed (in fact, La- 
grangian) coordinate in Eq. (29) radically differs from the 
form suggested in Refs. 5,7 on phenomenological grounds. 

To illustrate Eq. (29), which is one of the main results 
of the paper, we present an explicit expression for the 
x-only stress tensor 




where G'^(p, n) = — fkfk+p the exchange pair cor- 
relation function in a homogeneous system with the den- 
sity n. For a 3D system with Coulomb interaction from 
Eq. (31) follows 



where F{y) = [1 - 3?y/2 + y^/2]9{l - y). It important to 
note that this equation (and more general Eq. (31)) is, 
in fact, the explicit small- gradient limit of the dynamic 
OEP functional. 



A construction of the explicit form of AP^^ beyond the 
x-only approximation requires a knowledge of the Lan- 
dau functional or, more precisely, the functional P^^\n-p\ 
of Eq. (11) for a homogeneous system. This is obviously 
in a spirit of LDA, which treats an inhomogeneous system 
as a locally homogeneous one. In principle, it should be 
possible to calculate Pixv[np] using the non-equilibrium 
Keldysh technique^^. As is evident from Eq. (11), the 
tensor P;j;y[np] is uniquely defined by the quasiparticle 
energies Ep and the total energy E. The both quantities 
are the functionals of the full Green's function G. Using 
the general reconstriiction formulas^^ one can relate G to 
the quasiparticle distribution rip and thus restore the re- 
quired functional. If the system evolution starts from the 
ground state, the derived dynamical local approximation 
actually requires the functional Pfiu[np\ only on a set 
of t- and x-independent functions rip which describe a 
Fermi-type occupation of a volume of the p-space that is 
enclosed in a second-order surface, as given by Eq. (22) 
for iV(p,x) = 9{kp{x) — p^). This restriction should 
essentially simplify a calculation of P^vlup], which still 
remains a demanding, though feasible, task. 



In conclusion, we derived the nonlinear non-adiabatic 

counterpart of LDA. We showed that the xc stress ten- 
sor (and consequently the xc vector potential) is a local 
functional of one vector and one tensor variables, which 
describe dynamical deformation of coordinate and mo- 
mentum space respectively. 
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